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I .    Introduction 

This  report  describes  an  analysis  of  a  model  of  the  active 
sonar  detection  process.   The  analysis  was  done  in  order  to 
estimate  the  effects  of  signal  excess  fluctuations  on  the 
probability  of  detecting  a  sonar  echo.   The  analysis  is  an 
extension  of  an  analysis  that  is  described  in  Reference  1. 

The  model  has  the  following  characteristics:   Each  return 
corresponding  to  a  resolution  cell  has  an  associated  signal 
excess.   A  return  is  conditionally  classified  as  an  echo  if  and 
only  if  the  signal  excess  associated  with  the  return  is  greater 
than  or  equal  to  zero  (a  success) .   A  return  is  unconditionally 
classified  as  an  echo  if  and  only  if  it  is  the   k*1*1   success  in  a 
sequence  of  no  greater  than   n   returns.   This  detection 
criterion  has  been  referred  to  as  a  k-out-of-n  criterion. 

For  a  resolution  cell  that  does  not  contain  a  target,  the 
signal  excess  is  negative  infinity.   For  a  resolution  cell  that 
does  contain  a  target,  the  signal  excess  is  a  random  quantity 
that  is  determined  by  a  stochastic  process.   Three  stochastic 
processes  were  considered  in  the  analysis.   The  first  is  the 
gauss-markov  process,  the  second  is  the  lambda-sigma  jump  process 
and  the  third  is  a  modification  of  this  process.   In  the  analysis 
that  is  described  in  Reference  1,  only  the  lambda-sigma  jump 
process  and  its  modification  were  considered.   A  discussion  of 
the  gauss-markov  process  which  supports  its  inclusion  in  the 
analysis  is  contained  in  Reference  2  and  an  example  of  its  use  in 
an  active  sonar  detection  model  is  contained  in  Reference  3. 


The  active  sonar  detection  models  that  are  described  in 
Reference  3,  Reference  4  and  Reference  5  are  all  based  on  a 
k-out-of-n   criterion.    In  Reference  5,  it  is  suggested  that  k 
equal  to  three  and   n   equal  to  five  give  an  adequate  description 
of  the  active  sonar  detection  process  in  which  an  operator 
determines  the  classification  of  the  sonar  returns.   In  Reference 
6,  arguments  based  on  operational  data  are  given  both  for 
determining  sonar  return  classification  in  terms  of  signal  excess 
and  for  a   k-out-of-n   criterion. 

Since  each  return  corresponding  to  a  resolution  cell  has  an 
associated  signal  excess,  the  model  can  accommodate  false 
contacts  based  on  signal  inputs  but  not  false  contacts  based  on 
noise  alone.   However,  if  in  the  model  the  probability  of 
conditionally  classifying  noise  as  an  echo  were  .001,  then  with  a 
3-out-of-5   criterion  the  chance  of  a  false  contact  based  on 
noise  alone  would  be  less  than  one  in  ten  million. 


II.   The  Active  Sonar  Detection  Model 

In  the  model  of  the  active  sonar  detection  process,  the 
classification  of  an  echo  is  determined  by  its  signal  excess.   A 
sonar  return  is  conditionally  classified  as  an  echo  if  and  only 
if  the  signal  excess  is  greater  than  or  equal  to  zero.   The 
signal  excess  is  the  difference  between  the  echo  signal-to-noise 
ratio  in  decibels  and  the  sonar  receiver  recognition 
differential,  both  of  which  are  independent  random  variables  in 
the  model.   With   XgE(t)   the  signal  excess  for  an  echo  generated 
at  a  time  indicated  by  the  index  t,  the  signal  excess  is: 

(1)  XSE(t)  =  10  log(S/N)  -  XRD(t) 

where   10  log(S/N)   is  the  signal-to-noise  ratio  in  decibels  for 
an  echo  generated  at  that  time  and  XRD(t)   is  the  recognition 
differential  for  the  echo.   The  random  signal-to-noise  ratio  in 
decibels  is  determined  by  the  following  sonar  equation: 

(2)  10  log(S/N)  =  XSL(t)  -  2  XTL(t)  +  XTS(t)  -  [XNL(t)  -  XDI(t)] 
where,   XslC^)   is  the  sonar  source  level,   XTL(t)   is  the  one- 
way transmission  loss  between  the  target  and  the  sonar,   X-jis(t) 
is  the  target  strength,   XNL(t)   is  the  noise  level  and   XDj(t) 
is  the  sonar  directivity  index.   The  random  variable   XNL(t) 
accounts  for  ambient  noise,  self  noise  and  reverberation  and 
represents  a  power  sum  of  these  quantities. 

From  Equations  1  and  2,  the  mean  value  of  the  signal  excess 
is  given  by: 

(3)  SE(t)  =  SL(t)  -  2  TL(t)  +  TS(t)  -  [NL(t)  -  DI(t)]  -  RD(t) 
where  each  term  represents  the  mean  value  of  its  corresponding 


random  variable.   In  the  model,  the  mean  values  are  the  values 
that  would  be  used  in  the  sonar  equation  that  Equation  3  usually 
represents.   From  Equation  1  and  Equation  3  the  signal  excess  can 
be  expressed  as  follows: 
(4)    XSE(t)  =  SE(t)  +  X(t) 

where   X(t)   is  a  random  variable  whose  mean  is  zero  and  whose 
variance  is  equal  to  the  variance  of   XsE(t).   Equation  4  implies 
that  for  each  index  t  that  is  involved  in  describing  an 
encounter  there  is  a  random  variable   X(t)   that  determines  the 
fluctuations  in  the  signal  excess  at  the  time  corresponding  to 
the  index  t. 

The  three  stochastic  processes  that  were  used  in  the 
analysis  to  determine  the  random  variables  represented  by  X(t) 
are  described  in  Section  III.   Unfortunately,  it  appears  that  the 
following  quotation  from  Reference  7  has  some  relevance  for  each 
of  the  processes:   "As  discussed  in  Volume  I,  it  is  common 
practice  in  performance  analyses  to  model  signal  or  noise  or 
signal  excess  fluctuations  as  stochastic  processes.   The  choice 
of  a  specific  process  may  be  based  on  experience  but  usually 
tends  toward  mathematical  convenience." 


III.  Two  Simulation  Programs 

Two  simulation  programs  for  an  IBM-PC  were  used  in  the 
analysis  of  the  active  sonar  detection  model.   The  two  programs 
differ  in  the  stochastic  process  that  determines  the  random 
components  of  the  signal  excess.   In  the  first  program,  the 
stochastic  process  is  either  a  lambda-sigma  jump  process  or  a 
modification  of  a  lambda-sigma  jump  process.   In  the  second 
program,  the  stochastic  process  is  a  gauss-markov  process.   The 
programs  which  are  written  in  BASICA  are  described  in  Appendix  1 
and  are  listed  in  Appendix  2. 

The  programs  simulate  encounters  between  a  submarine  and  an 
active  sonar.   In  an  encounter,  the  sonar  source  level,  the 
transmission  loss,  the  submarine  target  strength,  the  noise  level 
and  the  sonar  directivity  index  are  independent  of  azimuth.   In 
addition,  the  course,  speed  and  depth  of  the  submarine  and  the 
depth  of  the  transmitting  and  receiving  hydrophones  of  the  sonar 
are  constant  throughout  an  encounter.   This  implies  that  the 
track  of  the  submarine  relative  to  the  sonar  is  a  straight  line 
and  that  the  encounter  can  be  determined  by  specifying  following 
quantities:  the  depth  of  the  submarine,  the  depth  of  the 
transmitting  hydrophone,  the  depth  of  the  receiving  hydrophone, 
the  relative  speed  of  the  submarine,  the  horizontal  range  of  the 
submarine  at  its  closest  point  of  approach  (CPA) ,  the  time  at  the 
beginning  of  the  encounter  and  the  time  at  the  end  of  the 
encounter.   An  encounter  of  this  kind  that  begins  and  ends  at 
ranges  beyond  which  the  probability  of  detecting  the  submarine  is 


essentially  zero  is  called  a  complete  straight  line  encounter. 
For  a  complete  straight  line  encounter,  only  the  horizontal  range 
at  CPA  is  required  to  determine  the  encounter  geometry.   This 
range  is  called  the  lateral  range.   The  probability  that 
detection  occurs  during  a  complete  straight  line  encounter  can  be 
expressed  as  a  function  of  the  lateral  range.   The  function  is 
called  the  lateral  range  function  or  lateral  range  curve. 

The  lateral  range  functions  plots  that  are  in  Section  IV 
were  generated  with  the  two  simulation  programs.   The  mean 
transmission  loss  values  that  were  used  to  do  this  are  plotted  in 
Figure  1  of  Section  IV.   The  values  correspond  to  horizontal 
ranges  from  2000  meters  to  100,000  meters  in  2000  meter 
increments.   For  a  target  range  that  is  not  an  integer  multiple 
of  2000  meters,  the  mean  transmission  loss  value  that  is  used  in 
the  programs  is  the  value  corresponding  to  the  greatest  integer 
multiple  that  is  less  than  the  target  range.   Relative  to  the 
programs,  the  transmission  loss  is  determined  by  a  step  function 
with  a  step  every  2000  meters. 

Because  of  the  step  function  nature  of  the  transmission 
loss,  the  relative  track  of  the  target  in  an  encounter  is  divided 
into  sectors  such  that  when  the  target  is  within  a  sector  the 
transmission  loss  is  constant.   The  simulation  programs  generate 
estimates  of  both  the  conditional  and  the  unconditional 
probability  that  detection  occurs  on  a  sector. 

Clearly,  the  encounters  represented  by  the  simulation 
programs  are  abstractions.   In  particular,  except  for  the  mean 


transmission  loss   TL(t) ,   all  terms  on  the  right  side  of 
Equation  2  are  constant  during  an  encounter.   Changes  due  to 
orientation,  location  and  relative  motion  are  ignored.   For 
example,  in  addition  to  ignoring  the  change  in  the  noise  field 
with  azimuth,  the  change  in  the  noise  field  due  to  the  doppler 
gain  that  results  from  the  radial  motion  of  the  sonar  and  the 
submarine  is  also  ignored.   (More  detailed  models  are  described 
in  Reference  3,  Reference  4  and  Reference  5.)   However,  the 
simulation  programs  were  assumed  to  be  adequate  to  satisfy  the 
purpose  of  the  analysis. 

Because  of  the  conditions  described  above,  Equation  3  can  be 
written  as: 

(5)   SE(t)  =  SL  -  2  TL[r(t)]  +  TS  -  (NL  -  DI)  -  RD 
where  the  transmission  loss  at  a  time  corresponding  to  the  index 
t   is  determined  by  the  sonar  transmitting  hydrophone  depth,  the 
sonar  receiving  hydrophone  depth,  the  submarine  depth  and   r(t), 
the  submarine  horizontal  range  at  that  time. 

Both  the  lambda-sigma  jump  process  and  the  gauss-markov 
process  have  been  used  to  predict  operational  performance  of 
sonar  systems.   This  is  not  the  case  for  the  modified  lambda- 
sigma  jump  process.   The  modification  was  introduced  in  order  to 
deal  with  the  observations  that  are  reported  in  Reference  8.   A 
modification  based  on  a  log-normal  distribution  rather  than  a 
shifted  rayleigh  distribution  would  have  accomplished  this  also. 
The  choice  of  a  rayleigh  distribution  was  based  on  mathematical 
convenience. 


For  a  rayleigh  random  variable  Y,  the  density  function  of 
its  distribution  is:   fy(y)  =  2  *a*yexp  (-a»y2 )   where   y  >  0. 
With   (3  =  T(3/2)/  the  mean  of  the  distribution  is:   jiy  =  P/di 
and  the  variance  is:   ay  =  (1  -  (32)/<x. 

The  modified  lambda-sigma  jump  process  is  defined  as 
follows:   If  X(t)   is  a  random  variable  that  is  determined  by  a 
modified  lambda-sigma  jump  process,  then  X(t)  =  p-y  -  YN(t) 
where  N(t)   is  determined  by  a  Poisson  process  with  a  mean  rate 
lambda  and  the  random  variables   Y^,  Y2 ,  '•'   are  independent 
rayleigh  random  variables  each  with  parameter  a.   This  implies 
that   p.x  =  0   and  <?x  =  aY-   If  a  =  (1  -  P:)*(l/aJ),  then 
ax  =  a.   Consequently,  if  the  random  component  of  the  signal 
excess  is  determined  by  a  modified  lambda-sigma  jump  process, 
then  the  standard  deviation  of  the  signal  excess  can  be  made 
equal  to  sigma  by  requiring  that   a  =  (1  -  (32)«(l/a2). 

The  correlation  coefficient  between  two  random  variables  is 
a  measure  of  their  dependence.   For  the  random  variables  X(t^) 
and  X(tj)   which  represent  the  random  component  of  the  signal 
excess  at  times   tj_   and  tj   and  for  tj  >  t^,   the  correlation 
coefficient  has  the  form  exp[-(tj  -  tjJ/T]   for  the  gauss-markov 
process,  the  lambda-sigma  jump  process  and  the  modified  lambda- 
sigma  jump  process.   For  each  process,  a  value  for  the  parameter 
t   and  a  value  for  the  parameter  a      are  sufficient  to  define  the 
process.   The  parameter   t   is  referred  to  as  the  relaxation 
time.   It  is  the  reciprocal  of  the  parameter  lambda  in  the 
lambda-sigma  jump  processes.   For  both  the  lambda-sigma  jump 
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jump  processes.   For  both  it  and  the  modified  lambda-sigma  jump 
process,  t   is  the  expected  time  between  jumps.   If   t  =  0   for 
any  one  of  the  processes,  the  random  components  of  the  signal 
excess  are  independent.   Otherwise,  the  random  components  are 
dependent.   If  (1/t)  =  0   for  any  of  the  processes,  the  random 
components  are  completely  dependent.   This  implies  that  the 
random  components  of  the  signal  excess  associated  with  an 
encounter  all  have  the  same  value.   Consequently,  if  the  value  of 
one  of  the  random  components  of  the  signal  excess  is  determined 
at  some  point  in  the  encounter,  the  value  of  all  of  the  other 
random  components  are  also  determined.   Reference  9  contains  a 
discussion  of  factors  affecting  relaxation  times. 

An  acceptable  choice  for   a,  the  standard  deviation  of  the 
signal  excess,  appears  to  be  more  easily  agreed  upon  than  an 
appropriate  choice  for   t,  the  relaxation  time.   However,  the 
following  quotation  from  Reference  10  is  relevant  both  to  these 
choices  and  to  the  choice  of  an  appropriate  stochastic  process: 
"As  of  this  date  [1974]  there  is  no  unanimity  of  opinion  as  to 
which  process  [lambda-sigma  jump  or  gauss-markov]  to  use  in  a 
detection  model  or  what  values  of  the  parameters  lambda  and  sigma 
to  use.   It  is  not  even  clear  (certainly  not  to  everyone)  that 
this  kind  of  detection  model  is  uniformly  valid.   Some  half- 
hearted attempts  at  'validating'  this  kind  of  model  have  been 
made  but  the  results  can  be  characterized  as  at  best  fragmentary 
and  inconclusive.   In  my  opinion,  what  is  needed  is  a  full-blown 
exercise  whose  sole  purpose  is  to  investigate  the  validity  of 


various  hypothesized  models  for  detection.   Past  attempts  at  this 
have  merely  been  'ex  post  facto'  exhumations  of  some  sketchy 
exercise  data." 
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IV.   Simulation  Results 

Encounter  simulation  data  are  presented  in  this  section. 
The  data  which  are  presented  as  lateral  range  function  plots  and 
as  range  sector  probability  of  detection  tables  are  from 
simulations  that  are  based  on  the  propagation  loss  values  that 
are  plotted  in  Figure  1,  the  encounter  values  that  are  listed  in 
Table  1  and  a   3-out-of-5   criterion. 

For  the  simulated  encounters,  the  lateral  range  function 
plots  in  Figure  2  through  Figure  6  indicate  the  effect  on  the 
encounter  probability  of  detection  of  the  transmission  loss 
function,  especially  the  increase  with  range  because  of  the 
convergence  zone  component.   But,  in  particular,  the  figures 
indicate  the  effect  on  the  encounter  'probability  of  detection  of 
a  gauss-markov  process  relative  to  that  a  lambda-sigma  jump 
process  for  various  relaxation  times.   Figure  7  is  representative 
of  the  effect  on  encounter  detection  probability  of  a  modified 
lambda-sigma  jump  process  relative  to  that  of  the  unmodified 
process.   For  additional  comparison  plots  for  the  same  simulated 
encounters,  see  Reference  1. 

Except  for  the  CPA  range  sector,  each  range  sector  is 
associated  with  two  track  segments  that  are  symmetric  about  the 
CPA.   The  CPA  range  sector  is  associated  with  a  single  track 
segment  that  is  bisected  by  the  CPA.   Table  2  and  Table  3  list 
the  probability  of  detection  estimates  on  a  track  segment  as  a 
function  of  the  track  segment's  distance  from  CPA  and  the  range 
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sector  associated  with  the  track  segment.   Table  4  and  Table  5 
do  this  for  the  conditional  probability  of  detection  estimates. 
The  values  in  Table  2  and  Table  3  are  estimates  of  the 
probability  that  detection  will  occur  on  a  track  segment.   The 
values  in  Table  4  and  in  Table  5  are  estimates  of  the  probability 
that  detection  will  occur  on  a  track  segment  given  detection  has 
not  already  occurred  on  an  earlier  track  segment. 

The  length  of  the  simulated  straight  line  encounters  that 
were  used  in  the  analysis  was  restricted  by  the  requirement  that 
the  range  from  the  sonar  at  the  end  points  of  an  encounter  be 
less  than  or  equal  to  100  kilometers.   Because  of  this 
restriction,  the  encounter  track  length  is  12  0  kilometers  for  an 
encounter  lateral  range  of  80  kilometers  and  the  encounter  track 
length  is  approximately  2  00  kilometers  for  an  encounter  lateral 
range  of  10  kilometers.   However,  for  the  simulated  acoustic 
conditions,  the  transmission  loss  for  a  range  of  100  kilometers 
is  94.2  dB.   If  transmission  loss  values  for  ranges  greater  than 
100  kilometers  were  constant  and  equal  to  94.2  dB  rather  than 
generally  increasing,  for  an  echo  from  a  target  that  was  at  a 
range  greater  than  100  kilometers,  the  probability  of 
conditionally  classifying  the  echo  as  an  echo  would  be  in  the 
neighborhood  of  .000015.   Thus,  increasing  the  length  of  the 
simulated  encounters  should  only  make  a  negligible  increase  in 
the  value  of  the  encounter  probability  of  detection  estimates. 

The  probability  estimates  generated  by  the  programs 
represent  samples  from  binomial  distributions  to  the  degree  of 
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randomness  of  the  program  simulations.   However,  the  sample  size 
of  the  estimates  is  not  constant.   In  a  simulated  encounter,  the 
encounter  is  terminated  when  a  detection  occurs.   Consequently, 
for  a  given  number  of  repetitions  of  an  encounter  simulation,  the 
sample  size  for  a  range  sector  depends  on  the  number  of 
detections  that  have  occurred  in  track  segments  prior  to  the 
track  segment  associated  with  the  range  sector.   Because  of  this, 
each  simulated  encounter  was  repeated  1000  times  in  order  to 
generate  the  lateral  range  function  plot  and  sector  probability 
of  detection  estimates.   This  number  of  samples  is  statistically 
sufficient  for  the  purpose  of  the  analysis  reported  here. 
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Figure  1.  A  plot  of  the  encounter  transmission  loss  values  in  dB 
against  their  corresponding  ranges  in  kilometers.  The  values  are 
listed  in  Appendix  4. 
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Recognition  Differential  20  dB 

Source  Level  215  dB 

Target  Strength  5  dB 

Noise  Level  65  dB 

Directivity  Index  20  dB 

Standard  Deviation:  a                                                                8  dB 

Ping  Cycle  Time  61  seconds 


Table  1.    The  parameter  values  for  the  simulated  encounters 
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Figure  2 .   Lateral  range  function  plots  for  a  gauss -markov 
process  and  a  lambda-sigma  jump  process  for  a  relaxation  time 
t   of  3  minutes. 
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Figure  3 .   Lateral  range  function  plots  for  a  gauss-markov 
process  and  a  lambda-sigma  jump  process  for  a  relaxation  time 
t   of  10  minutes. 
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Figure  4 .   Lateral  range  function  plots  for  a  gauss-markov 
process  and  a  lambda -sigma  jump  process  for  a  relaxation  time 
t   of  60  minutes. 
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Figure  5.   Lateral  range  function  plots  for  a  gauss-markov 
process  and  a  lambda -sigma  jump  process  for  a  relaxation  time 
t   of  zero. 
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Figure  6.   Lateral  range  function  plots  for  a  gauss-markov 
process  and  a  lambda -sigma  jump  process  for  a  relaxation  time 
t   of  infinity. 
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Figure  7.   Lateral  range  function  plots  for  a  lambda-sigma  jump 
process  and  a  modified  lambda-sigma  jump  process  for  a  relaxation 
time   t   of  10  minutes. 
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Table  2.   Range  sector  probability  of  detection  estimates  for  a 
gauss-markov  process,  a  lateral  range  of  10,000  meters,  a 
relaxation  time   t   of  3  minutes  and  a   3-out-of-5   criterion. 
Range  sector  and  distance  from  CPA  values  are  in  meters  and  zero 
probability  values  are  not  listed. 
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probability  range  sector  distance  from  CPA 

.01  76000  -74330 

.132  74000  -72312 

.089  72000  -70292 

.062  70000  -68272 

.029  68000  -66250 

.045  66000  -64226 

.007  64000  -62201 

.003  62000  -60175 

.001  50000  -47969 

.001  48000  -45924 

.001  40000  -37697 

.002  38000  -35623 

.001  26000  -22913 

.007  12000  -4583 

.007  10000  0 

.001  12000  4583 

.001  32000  29343 

.001  38000  35623 

.001  46000  43875 

.002  50000  47969 

.001  52000  50010 

.001  58000  56116 

.002  60000  58146 

.007  62000  60175 

.040  64000  62201 

.030  66000  64226 

.046  68000  66250 

.060  70000  68272 

.054  72000  70292 

.009  74000  72312 
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Table  3.   Range  sector  probability  of  detection  estimates  for  a 
lambda-sigma  jump  process,  a  lateral  range  of  10,000  meters,  a 
relaxation  time   t   of  3  minutes  and  a   3-out-of-5   criterion. 
Range  sector  and  distance  from  CPA  values  are  in  meters  and  zero 
probability  values  are  not  listed. 
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probability  range  sector  distance  from  CPA 

.01  76000  -74330 

.147  74000  -72312 

.090  72000  -70292 

.061  70000  -68272 

.046  68000  -66250 

.033  66000  -64226 

.017  64000  -62201 

.005  62000  -60175 

.001  60000  -58146 

.003  52000  -50010 

.002  50000  -47969 

.001  48000  -45924 

.003  46000  -43875 

.002  42000  -39762 

.001  40000  -37697 

.001  38000  -35623 

.001  28000  -25080 

.001  26000  -22913 

.001  24000  -20712 

.003  20000  -16155 

.002  18000  -13748 

.012  12000  -4583 

.008  10000  0 

.001  14000  8307 

.001  16000  11180 

.003  20000  16155 

.001  22000  18466 

.001  26000  22913 

.003  28000  25080 

.002  30000  27221 

.001  32000  29343 

.001  36000  35541 

.001  40000  37696 

.001  44000  41821 

.003  46000  43875 

.004  48000  45924 

.002  50000  47969 

.001  56000  54083 

.005  58000  56116 

.005  60000  58146 

.010  62000  60175 

.040  64000  62201 

.029  66000  64226 

.039  68000  66250 

.043  70000  68272 

.046  72000  70292 

.005  74000  72312 
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Table  4.   Range  sector  conditional  probability  of  detection 
estimates  for  a  gauss-markov  process,  a  lateral  range  of  10,000 
meters,  a  relaxation  time  of  3  minutes  and  a   3-out-of-5 
criterion.   Range  sector  and  distance  from  CPA  values  are  in 
meters  and  zero  probability  values  are  not  listed. 
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probability  range  sector  distance  from  CPA 

.01  76000  -74330 

.133  74000  -72312 

.104  72000  -70292 

.081  70000  -68272 

.041  68000  -66250 

.066  66000  -64226 

.011  64000  -62201 

.005  62000  -60175 

.002  50000  -47969 

.002  48000  -45924 

.002  40000  -37697 

.003  38000  -35623 

.002  26000  -22913 

.011  12000  -4583 

.011  10000  0 

.002  12000  4583 

.002  32000  29343 

.002  38000  35623 

.002  46000  43875 

.003  50000  47969 

.002  52000  50010 

.002  58000  56116 

.003  60000  58146 

.012  62000  60175 

.068  64000  62201 

.055  66000  64226 

.089  68000  66250 

.128  70000  68272 

.132  72000  70292 

.025  74000  72312 
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Table  5.   Range  sector  conditional  probability  of  detection 
estimates  for  a  lambda-sigma  jump  process,  a  lateral  range  of 
10,000  meters,  a  relaxation  time   t   of  3  minutes  and  a 
3-out-of-5   criterion.   Range  sector  and  distance  from  CPA  values 
are  in  meters  and  zero  probability  values  are  not  listed. 
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probability  range  sector  distance  from  CPA 

.01  76000  -74330 

.148  74000  -72312 

.107  72000  -70292 

.081  70000  -68272 

.066  68000  -66250 

.051  66000  -64226 

.027  64000  -62201 

.008  62000  -60175 

.002  60000  -58146 

.005  52000  -50010 

.003  50000  -47969 

.002  48000  -45924 

.005  46000  -43875 

.003  42000  -39762 

.002  40000  -37697 

.002  38000  -35623 

.002  28000  -25080 

.002  26000  -22913 

.002  24000  -20712 

.003  20000  -16155 

.005  18000  -13748 

.021  12000  -4583 

.014  10000  0 

.002  14000  8307 

.002  16000  11180 

.005  20000  16155 

.002  22000  18466 

.002  26000  22913 

.006  28000  25080 

.004  30000  27221 

.002  32000  29343 

.002  36000  35541 

.002  40000  37696 

.002  44000  41821 

.006  46000  43875 

.008  48000  45924 

.004  50000  47969 

.002  56000  54083 

.010  58000  56116 

.010  60000  58146 

.020  62000  60175 

.079  64000  62201 

.063  66000  64226 

.090  68000  66250 

.109  70000  68272 

.131  72000  70292 

.016  74000  72312 
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V.    Conclusions 

For  the  simulated  encounters,  the  lateral  range  function 
plots  in  Section  IV  suggest  that  the  probability  of  detection 
during  an  encounter  is  a  maximum  for  a  relaxation  time  in  the 
neighborhood  of  3  minutes  for  the  gauss-markov  process  as  well  as 
for  the  lambda-sigma  jump  process.   A  heuristic  argument  for  this 
apparent  relation  between  probability  of  detection,  relaxation 
time  and  detection  criterion  for  the  lambda-sigma  jump  process  is 
given  in  Appendix  3 . 

For  the  simulated  encounters  and  a  relaxation  time  of  3 
minutes,  the  lateral  range  function  plots  in  Section  IV  also 
suggest  that  the  probability  of  detection  during  an  encounter  for 
the  gauss-markov  process  and  for  lambda-sigma  jump  process  are 
essentially  equivalent.   This  apparent  equivalence  does  not  hold 
generally.   For  example,  for  relaxation  times  of  10  minutes  and 
60  minutes,  the  probability  of  detection  during  an  encounter  for 
the  gauss-markov  process  is  greater  than  that  for  the  lambda- 
sigma  jump  process.   However,  for  both  processes,  the  lateral 
range  function  plots  for  a  relaxation  time  of  zero  correspond  to 
encounters  in  which  the  random  components  of  the  signal  excess 
are  independent  gaussian  random  variables.   And,  for  both 
processes,  the  lateral  range  function  plots  for  a  relaxation  time 
of  infinity  correspond  to  encounters  in  which  the  random 
components  of  the  signal  excess  are  dependent  gaussian  random 
variables  with  correlation  coefficients  of  one.   Consequently,  in 
the  limit  for  a  relaxation  time  of  zero  or  infinity,  the  two 
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processes  are  equivalent.   Figure  5  shows  their  lateral  range 
plots  for  a  relaxation  time  of  zero  and  Figure  6  shows  their 
lateral  range  plots  for  a  relaxation  time  of  infinity. 

In  Table  6  on  Page  33,  for  the  nominal  encounter  ping  cycle 
time  of  60  seconds  and  a  lambda-sigma  jump  process,  conditional 
probability  of  detection  estimates  are  listed  for  a   1-out-of-l 
detection  criterion  and  a   3-out-of-5   detection  criterion.   (In 
a  1-out-of-l  detection  criterion,  a  return  is  unconditionally 
classified  as  an  echo  if  and  only  if  the  signal  excess  associated 
with  the  return  is  greater  than  or  equal  to  zero.) 

Also,  in  Table  6,  conditional  probabilities  of  detection  are 
listed  for  an  encounter  with  a  single  ping  per  sector,  1-out-of-l 
detection  criterion.   For  an  encounter  with  a  single  ping  per 
sector,  1-out-of-l  detection  criterion,   $(SE/a)  is  the 
conditional  probability  of  detection  for  a  range  sector  where   $ 
indicates  the  standard  normal  (gaussian)  cumulative  distribution 
function.   Conditional  probabilities  of  detection  determined  for 
a  range  sector  by   $(SE/a)   have  been  used  to  construct 
probability  of  detection  maps.   Table  6  suggests  that  for  a 
lambda-sigma  jump  process,  the  conditional  probabilities  of 
detection  for  a   3-out-of-5   detection  criterion  may  not  be 
significantly  different  from  the  values  determined  by   $(SE/a). 
However,  in  Appendix  3,  a  heuristic  argument  indicates  that  for  a 
lambda-sigma  jump  process,  a   3-out-of-5   encounter  conditional 
probability  of  detection  estimate  is  of  the  order  of 
$ (SE/a) • (t/T) *exp(-s/T)   where   s   is  the  time  required  for  3 
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consecutive  successes  to  occur  and  for  this  case   t   can  be 
considered  to  be  the  time  in  a  range  sector.   For  the  encounter 
simulation  that  generated  the  data  in  Table  6,  both  s  and   t 
are  equal  to  3  minutes  and,  for  the  range  sectors  with  reasonable 
sample  sizes,  t/T   is  approximately  equal  to  2.   This  implies 
that  the   3-out-of-5   values  should  be  of  the  order  of   $(SE/a) 
which  is  consistent  with  Table  6.   For  an  encounter  whose  lateral 
range  is  10  kilometers,  the  length  of  the  range  sector  track 
segments  go  from  zero  for  a  range  of  100  kilometers  to  13266 
meters  for  the  CPA  range  of  10  kilometers.   For  the  nominal  ping 
cycle  time  of  60  seconds,  the  number  of  pings  in  the  CPA  range 
sector  track  segment  is  approximately  42.   A  table  similar  to 
Table  6  is  generated  for  the  gauss-markov  process  for  the  same 
encounter  conditions. 

Table  6  also  suggests  that,  for  a  three  minute  relaxation 
time,  the   1-out-of-l  criterion  is   superior  to  the   3-out-of-5 
criterion  because  of  the  increased  probability  of  detection  and 
the  minimum  decision  delay.   However,  for  an  operational  sonar 
system,  the  false  alarm  probability  would  be  significantly  lower 
for  the   3-out-of-5   criterion  given  typical  single  return  false 
alarm  probabilities. 
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single  ping 

.000 
.017 
.130 
.130 
.077 
.070 
.067 
.032 
.008 
.002 
.001 
.001 
.002 
.004 
.005 
.005 
.003 
.002 
.002 
.002 
.002 
.002 
.002 
.001 
.002 
.002 
.002 
.001 
.001 
.001 
.002 
.002 
.002 
.007 


1-out-of-l 

3-out-of-5 

.047 

.015 

.323 

.127 

.192 

.107 

.180 

.079 

.119 

.055 

.149 

.064 

.028 

.022 

.026 

.010 

.003 

.002 

.007 

.002 

.007 

.007 

.002 

.017 

.005 

.003 

.005 

.017 

.003 

.003 

007 


007 
Oil 
003 

007 
Oil 
Oil 

004 

007 
007 
008 
111 


.002 

.002 
.002 


.002 


.005 


.005 
.022 


range  sector 

76000 

74000 

72000 

70000 

68000 

66000 

64000 

62000 

60000 

58000 

56000 

54000 

52000 

50000 

48000 

46000 

44000 

42000 

40000 

38000 

36000  . 

34000 

32000 

30000 

28000 

26000 

24000 

22000 

20000 

18000 

16000 

14000 

12000 

10000 


Table  6.   Range  sector  conditional  probabilities  of  detection  for 
a  single  ping  per  sector,   1-out-of-l   criterion  and  conditional 
probability  of  detection  estimates  for  a  lambda-sigma  jump 
process  and  a   1-out-of-l   criterion,  a   3-out-of-5   criterion 
and  a  10,000  meter  lateral  range.   Range  sector  values  are  in 
meters  and  zero  probability  values  are  not  listed. 
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Appendix  1.  Program  Descriptions 

The  programs  that  were  used  to  compute  encounter  detection 
probability  estimates  are  listed  in  Appendix  2.   The  programs  are 
written  in  BASICA  and  were  run  on  an  IBM-PC  using  QuickBasic 
which  is  licensed  by  Microsoft  Corporation. 

When  a  program  is  run,  an  initial  prompt  provides  a  choice 
of  either  generating  a  file  containing  encounter  probability 
estimates  or  printing  the  data  contained  in  a  previously 
generated  file.   When  the  program  is  used  to  generate  encounter 
probability  estimates,  additional  prompts  request  parameter 
values  that  are  required  to  define  an  encounter. 
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Appendix  2 .  Program  Listings 


file  (g/p)";N$ 

230 

8  0  ELSE  GOTO  50 


10  REM  active  gauss-markov  lateral  range  function/range  sector 

probabilities 

20  P$="RSGMA.BAS" 

30  DIM  TL(50) ,PP(20) ,N0(20) ,LR(20) ,CTA(20) ,M1(150) ,P1(150,20) , 

P2 (150,2  0) ,R1(150,2  0) ,¥1(150,20) 

40  PI=4*ATN(1) 

50  INPUT  "generate  file  or  print 

60  IF  N$="G"  OR  N$="g"  THEN  GOTO 

70  IF  N$="P"  OR  N$="p"  THEN  GOTO 

80  PRINT 

90  ON  ERROR  GOTO  100:  GOTO  110 

100  RESUME  90 

110  INPUT  "file  name";A$ 

120  OPEN  "I",#1,A$ 

13  0  INPUT# 1 , RD , SL , TS , NL, DI , TAU , SIG , W , LRM , ST , CT , N3 , NO , N8 , M8 , B$ , 

P$,L 

140  FOR  J=l  TO  L 

150  INPUTH,  PP(J)  ,LR(J)  ,M1(J)  ,N0(J)  ,CTA(J) 

160  FOR  1=0  TO  Ml (J) 

170  INPUT#1,  P1(I,J) ,R1(I,J) ,Y1(I,J) 

180  NEXT  I 

190  NEXT  J 

2  00  CLOSE 

210  ON  ERROR  GOTO  0 

220  GOTO  1440 

230  PRINT 

240  INPUT  "maximum  transmission  loss  range  (50000/100000) " ;MTLR 

250  IF  MTLR=50000!  THEN 

ELSE  GOTO  23  0 

260  INPUT  "transmission  loss  data  entry 

270  IF  A$="D"  OR  A$="d"  THEN  GOTO  400 

280  IF  A$="K"  OR  A$="k"  THEN  GOTO  290  ELSE  GOTO  260 

290  FOR  1=1  TO  50 

300  A$=STR$ (1*2000) 

310  PRINT  "transmission 

320  NEXT  I 

330  TL(0)=TL(1) 

340  INPUT  "transmission 

350  OPEN  "0",#1,B$ 

360  FOR  1=0  TO  50 

370  WRITE#1,  TL(I) 

380  NEXT  I 

390  CLOSE:  GOTO  490 

400  ON  ERROR  GOTO  410: 

410  RESUME  4  00 

420  INPUT  "transmission 

430  OPEN  "I",#1,B$ 

440  FOR  1=0  TO  50 

450  INPUT#1,  TL(I) 

460  NEXT  I 


RI=1000  ELSE  IF  MTLR=100000!  THEN  RI=200C 

(K=keyboard/D=disk) "  ;A$ 
3E  GOTO  260 

loss  at  "A$"  meters  (dB)":  INPUT  TL(I) 


loss  file  name";B$ 


GOTO  42  0 


loss  file  name";B$ 
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470  CLOSE 

480  ON  ERROR  GOTO  0 

490  INPUT  "recognition  differential  (dB) " ;RD 

500  INPUT  "source  level  (dB)";SL:  INPUT  "target  strength  (dB)";TS 

510  INPUT  "noise  level  (dB)";NL:  INPUT  "directivity  index  (dB)"; 

DI 

520  FOM=SL+TS-NL+DI-RD:  REM  figure  of  merit 

530  INPUT  "relaxation  time  tau  (minutes,  -1  for  infinity )" ;TAU 

540  INPUT  "sigma  (dB)";SIG 

550  INPUT  "relative  speed  (knots) ";W 

560  IF  W<=0  THEN  PRINT  "must  be  greater  than  0":  GOTO  550 

570  W=W*1852/60:  REM  relative  speed  in  meters  per  minute 

580  INPUT  "maximum  lateral  range  (meters) " ;LRM 

590  IF  LRM>.8*MTLR  THEN  PRINT  "maximum  is  " 

+STR$(.8*MTLR)+"  meters":  GOTO  580 

600  INPUT  "lateral  range  step  (meters) ";ST 

610  IF  ST>LRM  THEN  PRINT  "maximum  step  is  "+STR$(LRM):  GOTO  600 

620  L=INT(LRM/ST) 

630  IF  L>20  THEN  PRINT  "minimum  step  is  "+STR$ (LRM/20) :  GOTO  600 

640  INPUT  "ping  cycle  time  (seconds) " ;CT:  CTM=300 

650  IF  CT>CTM  THEN  PRINT  "must  be  less  than  "+STR$(CTM)+ 

"  seconds":  GOTO  640 

660  CT=CT/60:  REM  ping  cycle  time  in  minutes 

670  N5=INT(SQR(MTLR*MTLR-LRM*LRM)/W/CT) :  REM  number  of  pings  to 

CPA 

680  INPUT  "pings  in  detection  criterion  window" ;N8 

690  IF  N8>=2*N5  THEN  PRINT  "must  be  less  than  "+STR$ (2*N5) : 

GOTO  68  0 

700  IF  N8<1  THEN  PRINT  "minimum  is  1  :  GOTO  480 

710  INPUT  "echoes  required  for  detection" ;M8 

720  IF  M8<1  THEN  PRINT  "minimum  is  1":  GOTO  710 

730  IF  M8>N8  THEN  PRINT  "maximum  is  "+STR$(N8):  GOTO  710 

740  INPUT  "repetitions" ;N3 

750  IF  N3<=0  THEN  PRINT  "must  be  greater  than  0":  GOTO  740 

760  INPUT  "data  file  name";A$ 

770  LR=0 

780  FOR  J=l  TO  L 

790  LR=LR+ST:  LR(J)=LR:  REM  lateral  range 

800  Y0=SQR(MTLR*MTLR-LR*LR) :  REM  distance  from  CPA  at  encounter 

start 

810  TE=Y0/W:  REM  TE=time  to  CPA  in  minutes 

820  N0=INT(TE/CT) :  IF  N0=0  THEN  N0=1:  REM  adjusted  number  of 

pings  to  CPA 

83  0  CTA=TE/N0:  REM  adjusted  ping  cycle  time  in  minutes 

840  DIM  SE(2*N0) ,P(2*N0) ,Y(2*N0) /R(2*N0) ,N(2*N0) 

850  TK=-TE 

860  FOR  K=0  TO  NO 

870  R1=SQR(LR*LR+W*W*TK*TK) :  IN=INT(R1/RI)  REM  transmission  loss 

range  index 

880  SE=FOM-2*TL(IN) 

890  SE(K)=SE:  TK=TK+CTA:  R(K)=IN*RI  REM  transmission  loss  range 

for  FOM 
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900  NEXT  K 

910  FOR  K=N0+1  TO  2*N0 

920  SE(K)=SE(2*N0-K) :  R(K) =R(2*N0-K) 

930  NEXT  K 

940  P=0 

950  FOR  1=1  TO  N3 

960  DIM  M(N8-1) 

970  GOSUB  2120 

980  X=SIG*Q 

990  IF  TAU  =-1  GOTO  1020 

1000  IF  TAU=0  THEN  RHO=0:  FAC=1:  GOTO  1020 

1010  RHO=EXP(-CTA/TAU) :  FAC=SQR(l-RHO*RHO) :  REM  gauss-markov 

parameters 

1020  FOR  K=0  TO  2*N0 

1030  IF  K=0  OR  TAU=-1  GOTO  1060 

1040  GOSUB  2120 

1050  X=SIG*FAC*Q+RHO*X 

1060  XSE=SE(K)+X 

1070  IF  K<=N8-1  THEN  J1=K:  GOTO  1120 

1080  FOR  J1=0  TO  N8-2 

1090  M(J1)=M(J1+1) 

1100  NEXT  Jl 

1110  J1=N8-1 

1120  IF  XSE>=0  THEN  M(J1) =1 ' ELSE  M(J1)=0 

1130  S=0 

1140  FOR  J1=0  TO  N8-1 

1150  S=S+M(J1) 

1160  NEXT  Jl 

1170  IF  S>=M8  GOTO  1200 

1180  NEXT  K 

1190  GOTO  1210 

1200  P=P+1:  P(K)=P(K)+1 

1210  ERASE  M 

1220  NEXT  I 

1230  PP(J)=P/N3 

1240  P1(0,J)=P(0) :  R1(0,J)=R(0) :  Y1(0,J)=-Y0 

1250  FOR  1=1  TO  2*N0 

1260  IF  R(I)=R(I-1)  THEN  N(I)=N(I-1) :  PI (N(I) , J) =P(I) +P1 (N(I) , J) 

GOTO  1310 

1270  N(I)=N(I-1)+1:  PI (N(I) , J) =P(I) :  Rl (N(I) ,J)=R(I) 

1280  IF  LR>R1(N(I) ,J)  THEN  Yl (N (I) , J) =0 :  GOTO  1310 

1290  Y1(N(I) /J)=SQR(R1(N(I) ,J)*R1(N(I) ,J)-LR*LR) 

1300  IF  Y1(N(I) ,J)<ABS(Y1(N(I-1) ,J) )  THEN  Yl (N (I) , J) =-Yl (N (I) , J) 

1310  NEXT  I 

1320  M1(J)=N(2*N0) :  N0(J)=N0:  CTA(J) =CTA*60 

1330  ERASE  P,Y,R,N,SE 

1340  NEXT  J 

1350  OPEN  ,,0"/#l/A$ 

13  60  WRITE # 1 , RD , SL, TS , NL, DI , TAU , SIG , W , LRM , ST , CT , N3 , NO , N8 , M8 , 

B$,P$,L 

1370  FOR  J=l  TO  L 

1380  WRITE#1,  PP(J) ,LR(J) ,M1(J) ,N0(J) ,CTA(J) 
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1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 


(Y=yes/N=no) " ;D$ 


FOR  1=0  TO  Ml (J) 

WRITE* 1,  P1(I,J)  ,R1(I,J)  ,Y1(I,J) 

NEXT  I 

NEXT  J 

CLOSE 

PRINT 

INPUT  "print  encounter  parameters 

IF  D$="Y"  OR  D$="y"  GOTO  1480 

IF  D$="N"  OR  D$="n"  GOTO  1660  ELSE  GOTO  1450 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

PRINT 

INPUT  "print 

IF  D$="Y"  OR 

IF  D$="N"  OR 

LPRINT:  LPRINT: 

FOR  J=l  TO  L 

LPRINT :  LPRINT 

LPRINT  "lateral  range 

"probability  of  detection 
"adjusted  ping  cycle  time  (seconds) 


'program  file  name 

•transmission  loss  file  name 

'recognition  differential  (dB) 

'source  level  (dB) 

'target  strength  (dB) 

'noise  level  (dB) 

'directivity  index  (dB) 

'relaxation  time  tau  (minutes) 

•sigma  (dB) 

'relative  speed  (knots) 

'maximum  lateral  range  (meters) 

'lateral  range  step  (meters) 

'ping  cycle  time  (seconds) 

'pings  in  detection  criterion  window 

'echoes  required  for  detection 

'repetitions 

'data  file  name 


,P$ 
,B$ 
,RD 
,SL 

,TS 

,NL 

,DI 

,TAU 

,SIG 

,W*60/1852 

,LRM 

,ST 

/CT*60 

,N8 

,M8 

,N3 

,A$ 


encounter  probabilities  (Y=yes/N=no) " ;D$ 

D$="y"  GOTO  1700 

D$="n"  GOTO  1780  ELSE  GOTO  1670 

LPRINT:  LPRINT  A$+"  encounter  probabilities" 


»,LR(J) 
»,PP(J) 
",CTA(J) 


(Y=yes/N=no) ";D$ 


LPRINT 

LPRINT 

LPRINT 

NEXT  J 

PRINT 

INPUT  "print  range  sector  probabilities 

IF  D$="Y"  OR  D$="y"  GOTO  1820 

IF  D$="N"  OR  D$="n"  GOTO  1930  ELSE  GOTO  1790 

LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+"  range  sector 
probabilities" 
1830  FOR  J=l  TO  L 
1840  LPRINT:  LPRINT 

1850  LPRINT  "lateral  range  "/LR(J) 

1860  LPRINT  "probability  of  detection  ",PP(J) 

1870  LPRINT  "adjusted  ping  cycle  time  (seconds)     ",CTA(J) 
1880  LPRINT 
1890  FOR  1=0  TO  Ml (J) 
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1900  LPRINT  "P("I")  =  "P1(I,J)/N3  TAB(30)  "R("I")  =  "R1(I,J) 

TAB(60)  "Y("I")  =  "Y1(I,J) 

1910  NEXT  I 

1920  NEXT  J 

1930  PRINT 

1940  INPUT  "print  range  sector  conditional  probabilities 

(Y=yes/N=no) ";D$ 

1950  IF  D$="Y"  OR  D$="y"  GOTO  1970 

1960  IF  D$="N"  OR  D$="n"  GOTO  2110  ELSE  GOTO  1940 

1970  LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+"  range  sector  conditional 

probabilities" 

1980  FOR  J=l  TO  L 

1990  LPRINT:  LPRINT. 

2000  LPRINT  "lateral  range  "/LR(J) 

2010  LPRINT  "probability  of  detection  ",PP(J) 

2020  LPRINT  "adjusted  ping  cycle  time  (seconds)     ",CTA(J) 

2030  LPRINT 

2  040  D=N3 

2050  FOR  1=0  TO  Ml (J) 

2060  IF  D=0  THEN  P2(I/J)=0  ELSE  P2 (I , J) =P1 (I , J)/D 

2070  LPRINT  "P("I")  =  "P2(I,J)  TAB(30)  "R("I")  =  "R1(I,J)  TAB(60) 

"Y("I")  =  "Y1(I,J) 

2080  IF  D>0  THEN  D=D-P1(I,J) 

2090  NEXT  I 

2100  NEXT  J 

2110  END 

2120  RV=RND:  IF  RV=0  THEN  RV=2.8E-38 

2130  IF  RV=1  THEN  RV=l-2.8E-38 

2140  Q=SQR(-2*LOG(RV) ) *SIN (2*PI*RND) :  REM  standard  normal  random 

number 

2150  RETURN 
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10  REM  active  lambda-sigma  lateral  range  function/range  sector 

probabilities 

20  P$="RSLSA.BAS" 

30  DIM  TL(50) /PP(20) ,N0(20) ,LR(20) ,CTA(20) ,M1(150) ,P1(150,20) , 

P2(150/20) ,R1(150,20) ,¥1(150,20) 

40  PI=4*ATN(1):  RC1=. 88623:  RC2=. 21459 

50  INPUT  "generate  file  or  print  file  (g/p)";N$ 

60  IF  N$="G"  OR  N$="g"  THEN  GOTO  23  0 

70  IF  N$="P"  OR  N$="p"  THEN  GOTO  80  ELSE  GOTO  50 

80  PRINT 

90  ON  ERROR  GOTO  100:  GOTO  110 

100  RESUME  90 

110  INPUT  "file  name";A$ 

120  OPEN  "I",#l/A$ 

130  INPUT#1,  RD,SL,TS,NL,DI,TAU,SIG,W,LRM,ST,CT,N3,N0,N8,M8, 

B$,C$,P$,L 

140  FOR  J=l  TO  L 

150  INPUT* 1,  PP(J) /LR(J) ,M1(J) ,N0(J) ,CTA(J) 

160  FOR  1=0  TO  Ml (J) 

170  INPUT* 1,  P1(I,J) ,R1(I,J) ,Y1(I,J) 

180  NEXT  I 

190  NEXT  J 

200  CLOSE 

210  ON  ERROR  GOTO  0 

220  GOTO  1520 

23  0  PRINT 

240  INPUT  "maximum  transmission  loss  range  (50000/100000) " ;MTLR 

250  IF  MTLR=50000!  THEN  RI=1000  ELSE  IF  MTLR=100000!  THEN  RI=2000 

ELSE  GOTO  23  0 

260  INPUT  "transmission  loss 

270  IF  A$="D"  OR  A$="d"  GOTO 

280  IF  A$="K"  OR  A$="k"  GOTO 

290  FOR  1=1  TO  50 

300  A$=STR$ (1*2000) 

310  PRINT  "transmission  loss  at  "A$"  meters  (dB)":  INPUT  TL(I) 

320  NEXT  I 

330  TL(0)=TL(1) 

340  INPUT  "transmission  loss  file  name";B$ 

350  OPEN  "0",#1,  B$ 

360  FOR  1=0  TO  50 

370  WRITE#1/  TL(I) 

380  NEXT  I 

390  CLOSE:  GOTO  520 

400  ON  ERROR  GOTO  410:  GOTO  420 

410  RESUME  400 

420  INPUT  "transmission  loss  file  name";B$ 

430  OPEN  "I",#l,  B$ 

440  FOR  1=0  TO  50 

450  INPUT* 1,  TL(I) 

4  60  NEXT  I 

470  CLOSE 

480  ON  ERROR  GOTO  0 


data  entry  (K=keyboard/D=disk) " ;A$ 

400 

290  ELSE  GOTO  260 
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490  INPUT  "signal  excess  distribution  (G=gaussian/R=rayleigh) " ;C$ 

500  IF  C$="GM  OR  C$="g"  GOTO  520 

510  IF  C$="R"  OR  C$="r"  GOTO  520  ELSE  GOTO  490 

520  INPUT  "recognition  differential  (dB)";RD 

530  INPUT  "source  level  (dB)";SL:  INPUT  "target  strength  (dB)";TS 

540  INPUT  "noise  level  (dB)";NL:  INPUT  "directivity  index  (dB)"; 

DI 

550  FOM=SL+TS-NL+DI-RD:  REM  figure  of  merit 

560  INPUT  "relaxation  time  tau  (minutes,  -1  for  infinity) " ;TAU 

570  INPUT  "sigma  (dB)";SIG 

580  INPUT  "relative  speed  (knots) ";W 

590  IF  W<=0  THEN  PRINT  "must  be  greater  than  0":  GOTO  580 

600  W=W*1852/60:  REM  relative  speed  in  meters  per  minute 

610  INPUT  "maximum  lateral  range  (meters) " ;LRM 

62  0  IF  LRM>.8*MTLR  THEN  PRINT  "maximum  is  "+STR$ ( . 8*MTLR) + 

"  meters":  GOTO  610 

630  INPUT  "lateral  range  step  (meters) ";ST 

640  IF  ST>LRM  THEN  PRINT  "maximum  step  is  "+STR$(LRM):  GOTO  630 

650  L=INT(LRM/ST) 

660  IF  L>20  THEN  PRINT  "minimum  step  is  "+STR$ (LRM/2  0) +"  meters": 

GOTO  63  0 

670  INPUT  "ping  cycle  time  (seconds) " ;CT:  CTM=300 

680  IF  CT>CTM  THEN  PRINT  "must  be  less  than  "+STR$(CTM)+ 

"  seconds":  GOTO  670 

690  CT=CT/60:  REM  ping  cycle  time  in  minutes 

7  00  N5=INT(SQR(MTLR*MTLR-LRM*LRM)/W/CT) :  REM  number  of  pings  to 

CPA 

710  INPUT  "pings  in  detection  criterion  window" ;N8 

720  IF  N8>=2*N5  THEN  PRINT  "must  be  less  than  "+STR$ (2*N5) : 

GOTO  710 

730  IF  N8<1  THEN  PRINT  "minimum  is  1":  GOTO  670 

740  INPUT  "echoes  required  for  detection" ;M8 

750  IF  M8<1  THEN  PRINT  "minimum  is  1":  GOTO  740 

760  IF  M8>N8  THEN  PRINT  "maximum  is  "+STR$(N8):  GOTO  740 

770  INPUT  "repetitions" ;N3 

780  IF  N3<=0  THEN  PRINT  "must  be  greater  than  0":  GOTO  770 

790  INPUT  "data  file  name";A$ 

800  LR=0 

810  FOR  J=l  TO  L 

82  0  LR=LR+ST:  LR(J)=LR:  REM  lateral  range 

83  0  Y0=SQR(MTLR*MTLR-LR*LR) :  REM  distance  from  CPA  at  encounter 
start 

840  TE=Y0/W:  REM  TE=time  to  CPA  in  minutes 

850  N0=INT(TE/CT) :  IF  N0=0  THEN  N0=1:  REM  adjusted  number  of 

pings  to  CPA 

860  CTA=TE/N0:  REM  adjusted  ping  cycle  time  in  minutes 

870  DIM  SE(2*N0) ,P(2*N0) ,Y(2*N0) ,R(2*N0) ,N(2*N0) 

880  TK=-TE 

890  FOR  K=0  TO  NO 

900  R1=SQR(LR*LR+W*W*TK*TK) :  IN=INT (Rl/RI)  REM  TL  range  index 

910  SE=FOM-2*TL(IN) 
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920  SE(K)=SE:  TK=TK+CTA:  R(K)=IN*RI  REM  transmission  loss  range 

for  FOM 

930  NEXT  K 

940  FOR  K=N0+1  TO  2*N0 

950  SE(K)=SE(2*N0-K) :  R(K) =R(2*N0-K) 

960  NEXT  K 

970  P=0 

980  FOR  1=1  TO  N3 

990  DIM  M(N8-1) 

1000  IF  TAU=0  GOTO  1050 

1010  GOSUB  2210 

1020  IF  TAU=-1  GOTO  1050 

1030  GOSUB  2270 

1040  N=-TAU*LOG(l-RV)/CTA:  REM  number  of  pings  to  the  first  jump 

1050  FOR  K=0  TO  2*N0 

1060  IF  TAU=0  GOTO  1110 

1070  IF  K<N  OR  TAU=-1  GOTO  112  0 

1080  GOSUB  2270 

1090  DN=-TAU*LOG(l-RV)/CTA 

1100  N=N+DN:  REM  number  of  pings  to  the  next  jump 

1110  GOSUB  2210 

1120  XSE=SE(K)+Q 

1130  IF  K<=N8-1  THEN  J1=K:  GOTO  1180 

1140  FOR  J1=0  TO  N8-2 

1150  M(J1)=M(J1+1) 

1160  NEXT  Jl 

1170  J1=N8-1 

1180  IF  XSE>=0  THEN  M(J1)=1  ELSE  M(J1)=0 

1190  S=0 

1200  FOR  J1=0  TO  N8-1 

1210  S=S+M(J1) 

1220  NEXT  Jl 

1230  IF  S>=M8  THEN  GOTO  12  60 

1240  NEXT  K 

1250  GOTO  1270 

1260  P=P+1:  P(K)=P(K)+1 

1270  ERASE  M 

1280  NEXT  I 

1290  PP(J)=P/N3 

1300  N(0)=0:  P1(0,J)=P(0) :  Rl (0, J) =R(0) :  Y1(0/J)=-Y0 

1310  FOR  1=1  TO  2*N0 

1320  IF  R(I)=R(I-1)  THEN  N(I)=N(I-1) :  PI (N (I) , J) =P(I) +P1 (N(I) , J) 

GOTO  1370 

1330  N(I)=N(I-1)+1:  PI (N(I) , J) =P(I) :  Rl (N(I) , J) =R(I) 

1340  IF  LR>R1(N(I) , J)  THEN  Yl (N(I) , J) =0 :  GOTO  1370 

1350  Y1(N(I) /J)=SQR(R1(N(I) ,J)*R1(N(I) ,J)-LR*LR) 

1360  IF  Y1(N(I) ,J)<ABS(Y1(N(I-1) ,J) )  THEN  Yl (N(I) , J) =-Yl (N (I) , J) 

1370  NEXT  I 

1380  M1(J)=N(2*N0) :  N0(J)=N0:  CTA(J) =CTA*60 

1390  ERASE  P,Y,R,N,SE 

1400  NEXT  J 

1410  IF  C$="G"  OR  C$="g"  THEN  C$="gaussian" :  GOTO  1430 
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1420  C$="rayleigh" 

1430  OPEN  "0",#1,A$ 

1440  WRITE#1,  RD,SL/TS,NL/DI/TAU,SIG,W,LRM,ST,CT/N3,N0,N8, 

M8,B$,C$,P$,L 

1450  FOR  J=l  TO  L 

1460  WRITE#1,  PP(J) ,LR(J) ,M1(J) ,N0(J) ,CTA(J) 

1470  FOR  1=0  TO  Ml (J) 

1480  WRITE#1,  P1(I,J) ,R1(I,J) ,Y1(I,J) 

1490  NEXT  I 

NEXT  J 

CLOSE 

PRINT 

INPUT  "print  encounter  parameters  (Y=yes/N=no) 

IF  D$="Y"  OR  D$="y"  GOTO  1560 

IF  D$="N"  OR  D$="n"  GOTO  1750  ELSE  GOTO  1530 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

LPRINT 

PRINT 

INPUT  "print 

IF  D$="Y"  OR 

IF  D$="N"  OR 

LPRINT:  LPRINT: 

FOR  J=l  TO  L 

LPRINT:  LPRINT 

LPRINT  "lateral  range 

"probability  of  detection 
"adjusted  ping  cycle  time  (seconds) 


1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 


program  file  name 

transmission  loss  file  name 

signal  excess  distribution 

recognition  differential  (dB) 

source  level  (dB) 

target  strength  (dB) 

noise  level  (dB) 

directivity  index  (dB) 

relaxation  time  tau  (minutes) 

sigma  (dB) 

relative  speed  (knots) 

maximum  lateral  range  (meters) 

lateral  range  step  (meters) 

ping  cycle  time  (seconds) 

pings  in  detection  criterion  window 

echoes  required  for  detection 

repetitions 

data  file  name 


;D$ 


",P$ 
",B$ 
",C$ 
",RD 
",SL 
",TS 
",NL 


,DI 

,TAU 
,SIG 

,W*60/1852 
,LRM 
,ST 

,CT*60 
,N8 
,M8 
,N3 
",A$ 


ti 


ii 


encounter  probabilities  (Y=yes/N=no) " ;D$ 

D$="y"  GOTO  1790 

D$="n"  GOTO  1870  ELSE  GOTO  1770 

LPRINT:  LPRINT  A$+"  encounter  probabilities" 


LPRINT 

LPRINT 

LPRINT 

NEXT  J 

PRINT 

INPUT  "print  range  sector  probabilities 

IF  D$="Y"  OR  D$="y"  GOTO  1910 

IF  D$="N"  OR  D$="n"  GOTO  2020  ELSE  GOTO  1880 

LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+"  range  sector 


",LR(J) 
»,PP(J) 
",CTA(J) 


(Y=yes/N=no) ";D$ 


probabilities" 
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1920  FOR  J=l  TO  L 

1930  LPRINT:  LPRINT 

194  0  LPRINT  "lateral  range  " ,LR(J) 

1950  LPRINT  "probability  of  detection  ",PP(J) 

1960  LPRINT  "adjusted  ping  cycle  time  (seconds)     ",CTA(J) 

1970  LPRINT 

1980  FOR  1=0  TO  Ml (J) 

1990  LPRINT  "P("I")  =  "P1(I,J)/N3  TAB(30)  "R("I")  =  "R1(I,J) 

TAB(60)  "Y("I")  =  "Y1(I,J) 

2  000  NEXT  I 

2010  NEXT  J 

2020  PRINT 

2030  INPUT  "print  range  sector  conditional  probabilities 

(Y=yes/N=no) ";D$ 

2040  IF  D$="Y"  OR  D$="y"  GOTO  2060 

2050  IF  D$="N"  OR  D$="n"  GOTO  2200  ELSE  GOTO  2030 

2060  LPRINT:  LPRINT:  LPRINT:  LPRINT  A$+"  range  sector  conditional 

probabilities" 

2070  FOR  J=l  TO  L 

2080  LPRINT:  LPRINT 

2090  LPRINT  "lateral  range  ",LR(J) 

2100  LPRINT  "probability  of  detection  ",PP(J) 

2110  LPRINT  "adjusted  ping  cycle  time  (seconds)     ",CTA(J) 

2120  LPRINT 

2130  D=N3 

2140  FOR  1=0  TO  Ml (J) 

2150  IF  D=0  THEN  P2(I,J)=0  ELSE  P2 (I , J) =P1 (I, J)/D 

2160  LPRINT  "P("I")  =  "P2(I/J)  TAB(30)  "R("I")  =  "R1(I,J)  TAB(60) 

"Y("I")  =  "Y1(I,J) 

2170  IF  D>0  THEN  D=D-P1(I,J) 

2180  NEXT  I 

2190  NEXT  J 

2200  END 

2210  GOSUB  2270 

2220  IF  C$="G"  OR  C$="g"  GOTO  2250 

2230  Q=SIG*(RCl-SQR(-LOG(RV) ) )/SQR(RC2) :  REM  random  signal  excess 

part 

2240  RETURN 

2250  Q=SIG*SQR(-2*LOG(RV) ) *SIN(2*PI*RND) :  REM  random  signal 

excess  part 

2260  RETURN 

2270  RV=RND:  IF  RV=0  THEN  RV=2.8E-38 

2280  IF  RV=1  THEN  RV=l-2.8E-38 

2290  RETURN 
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Appendix  3 .  A  Heuristic  Argument 

The  following  heuristic  argument  applies  to  the  lambda-sigma 
jump  process.   Suppose  that   SE   is  constant  over  a  time  interval 
of  length  t.   Then,  the  occurrence  of  a  jump  during  the  interval 
that  is  sufficient  for  conditionally  classifying  an  echo  as  an 
echo  (a  success)  is  determined  by  a  Poisson  process  with  rate 
$(SE/<j)/t.   In  addition,  suppose  §(SE/a)   is  small  enough  so 
that  the  probability  that  one  success  will  occur  in  the  interval 
is  approximately   § (SE/a) • (t/T) .   (If  this  expression  is  equal  to 
.1,  the  probability  of  one  success  in  the  interval  is  .0905  and 
the  probability  of  two  or  more  successes  in  the  interval  is 
.0047.)   The  probability  that  the  time  to  the  first  jump 
following  a  success  will  be  greater  than  or  equal  to   s   is 
exp(-s/T).    This  implies  that  the  probability  of  a  success  in 
the  interval  after  which  the  random  component  of  the  signal 
excess  is  constant  for  a  time  at  least  equal  to   s   is 
approximately   $ (SE/a) • (t/T) •exp(-s/T) .   This  expression  is  a 
maximum  when   t  =  s.   Let   s   be  the  time  required  for  k 
consecutive  returns  and  let   t   be  a  time  greater  than  but  of  the 
order  of  the  time  required  for  n   consecutive  returns.   Then, 
given  that  the  signal  excess  at  the  beginning  of  a  time  interval 
of  length   t   is  less  than  zero,  the  above  suggests  that  for  a 
k-out-of-n   criterion  the  conditional  probability  of  detection 
during  the  interval   will  have  a  maximum  for   t   in  the 
neighborhood  of   s   if   $ (SE/a) • (t/T) »exp(-s/T)   is  small  and  s 
is  sufficiently  smaller  than   t.   This  suggests,   for   t   in  the 
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neighborhood  of   s,  that  the  range  sector  conditional 
probabilities  for  an  encounter  will  have  a  maximum  and 
consequently  the  detection  probability  for  the  encounter  will 
have  a  maximum  also. 

As   t   becomes  small,   $ (SE/a) • (t/T) »exp(-s/T)   will  become 
greater  than  one  and  the  approximation  described  above  clearly 
does  not  hold.   In  this  case,  given  the  signal  excess  at  the 
beginning  of  a  time  interval  of  length  t   is  less  than  zero,  the 
conditional  probability  of  detection  during  an  interval 
approaches  the  probability  of  3  or  more  successes  of  out  of  5 
independent  trials  where  the  probability  of  success  is   §(SE/a). 

As   t   becomes  large,  the  probability  of  the  random 
component  of  the  signal  excess  changing  during  an  encounter 
becomes  small  and  in  the  limit  the  random  component  of  the  signal 
excess  will  be  a  constant  value  during  an  encounter  and  the 
signal  excess  can  be  represented  by   SE  +  X  where  the  value  of 
X   is  determined  at  the  beginning  of  the  encounter.   In  this 
case,  if   SE  +  X   is  greater  than  zero  for  3  returns  out  of  5 
during  an  encounter,  detection  will  occur,  otherwise  it  will  not. 

For  the  encounters  that  were  considered,  for  both  small  and 
large  values  of   t,   the  conditional  probabilities  of  detection 
for  the  range  segments  of  an  encounter  appear  to  be  dominated  by 
those  for   t   in  the  neighborhood  of  3  minutes. 
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